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We present a measurement of the top quark mass with ft dilepton events produced in pp collisions 
at the Fermilab Tevatron (^=1-96 TeV) and collected by the CDF II detector. A sample of 328 
events with a charged electron or muon and an isolated track, corresponding to an integrated lumi- 
nosity of 2.9 fb _1 , are selected as it candidates. To account for the unconstrained event kinematics, 
we scan over the phase space of the azimuthal angles {4> V1 , <f> V2 ) of neutrinos and reconstruct the 
top quark mass for each <j> vi , 0„ 2 pair by minimizing a \ 2 function in the it dilepton hypothesis. 
We assign x 2 -dependent weights to the solutions in order to build a preferred mass for each event. 
Preferred mass distributions (templates) are built from simulated it and background events, and 
parameterized in order to provide continuous probability density functions. A likelihood fit to the 
mass distribution in data as a weighted sum of signal and background probability density functions 
gives a top quark mass of 165.5jl3;3(stat.)±3.1(syst.) GeV/c 2 . 



PACS numbers: 14.65.Ha, 13.85.Qk, 12.15.Ff 



'Deceased 



tWith visitors from a University of Massachusetts Amherst, 



4 



I. INTRODUCTION 

The standard model (SM) explains the non-zero weak 
boson masses by spontaneous breaking of the electroweak 
(EW) symmetry induced by the Higgs field Also, 
non-zero quark masses are generated by the coupling of 
the Higgs doublet with the fundamental fermions. How- 
ever, their values are not predicted since they are propor- 
tional to the unknown Yukawa couplings of each quark. 
The enormous top quark mass, which has a value com- 
parable to the EW scale, justifies the suspect that this 
quark may play a special role in the mechanism which 
breaks EW symmetry. In addition, because of its large 
mass, the top quark gives the largest contribution to loop 
corrections in the W propagator. Within the SM, the cor- 
relation between the top mass and the W mass induced 
by these corrections allows setting limits on the mass of 
the yet unobserved Higgs boson, and favor a relatively 
light Higgs. A more accurate measurement of the top 
quark mass will tighten the SM predicted region for the 
Higgs boson mass. 

According to the SM, at the Tevatron's 1.96 TeV en- 
ergy top quarks are dominantly produced in pairs, by 
qq annihilation in ~ 85% of the cases and by gluon fu- 
sion in the remaining ~ 15% [2J. Due to its extremely 
short lifetime, which in the SM is expected to be about 
1CP 25 s, the top quark decays before hadronizing in ~ 
100% of cases into a W boson and fr-quark Q. Subse- 
quently the W boson can either decay into quarks as a qq' 
pair or into a charged lepton-neutrino pair. This allows 
classifying the it candidate events into three final states: 
all-hadronic, lepton+jets, or dilepton, depending on the 
decay modes of the two W bosons in the event. The 
all-hadronic state, where both Ws decay hadronically 
(about 46% of it events) , is characterized by six or more 
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jets in the event. The lepton+jets final state contains one 
electron or muon (about 30% of it events), four or more 
jets, and one neutrino. Analyses dealing with the lep- 
ton+jets final state have provided the most precise top 
quark mass measurements, due to an optimal compro- 
mise between statistics and backgrounds. The dilepton 
final state, which is defined by the presence of two lep- 
tons (electrons or muons, about 5% of it events), two or 
more jets, and large missing transverse energy from the 
two neutrinos, is the cleanest one, but suffers from the 
poorest statistics. 

It is important to perform measurements using inde- 
pendent data samples in all final states in order to im- 
prove the precision on the top quark mass and to be able 
to cross-check the results. Once the channel-specific SM 
backgrounds have been removed, discrepancies in the re- 
sults across different samples could provide hints of new 
physics. The present analysis is performed in the dilep- 
ton final state by means of lepton+track ( "LTRK" ) top- 
pair selection. This selection is chosen to collect a large 
portion of events (about 45%) not involved in the other 
CDF high precision top mass analyses performed in the 
dilepton final state 0, [f| . 

The paper reports a measurement of the top quark 
mass with data collected by CDF II before spring 2008, 
corresponding to 2.9 fb _1 of integrated luminosity. We 
select it candidate events in dilepton channel by requir- 
ing a well identified electron or muon plus a second, more 
loosely defined lepton, which is an isolated track. The 
measurement of the top quark mass in this channel is 
particularly challenging because of the two neutrinos in 
the final state. The kinematics is under-constrained and 
therefore assumptions on some missing final state ob- 
servables are needed in order to reconstruct the event. 
In order to constrain the kinematics, we scan over the 
space of possibilities for the azimuthal angles of the two 
neutrinos, and reconstruct the top quark mass by mini- 
mizing a x 2 function using the it dilepton hypothesis. A 
weighted average over a grid of the azimuthal neutrino 
angles {<p Vl , 4> V2 ) returns a single top quark mass value 
per event. In this analysis the Breit-Wigner probability 
distribution function with a top quark mass-dependent 
decay width is applied in the kincmatical event recon- 
struction, which helps to decrease the statistical uncer- 
tainty by 20% compared to the method described in [||. 
The top quark mass distribution in the data is fitted to 
the parameterized signal and background templates, and 
the mass is extracted as the one corresponding to the 
best fit. 

Sections [TT] and IIIII describe the detector and the selec- 
tion of the data sample. Section IIVI gives an overview of 
the method used to reconstruct the events and to derive a 
single value of the top quark mass for each event. Section 
IVl defines the parameterization of signal and background 
mass distributions and the likelihood function used to fit 
the data to these distributions. Section [VTI describes the 
studies performed to calibrate the method, sections IVIII 
and I Villi present the results and the systematic uncer- 
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tainties, and section ITXl gives the conclusions. 



II. THE CDF II DETECTOR 

The Collider Detector at Fermilab was upgraded in 
the year 2000 (CDF II, Figure [J) in order to be able to 
handle the higher collision rate from the increased Teva- 
tron luminosity. CDF II is a cylindrically and forward- 
backward symmetric apparatus detecting the products of 
the pp collisions over almost the full solid angle. A cylin- 
drical (r, <fr, z) coordinate system is used to describe the 
detector geometry. The origin of the reference system 
is the geometric center of the detector, with the z axis 
pointing along the proton beam. The pseudorapidity 77 is 
defined by r\ = — ln(tan(#/2)), where 9 is the polar angle 
relative to the z axis. The detector elements which are 
most relevant for this analysis are described below. A 
more complete description of the detector can be found 
elsewhere QJ- 

The tracking system consists of an inner silicon sys- 
tem and an outer gas drift chamber, the Central Outer 
Tracker (COT). The entire tracker is enclosed in a su- 
perconducting solenoid which generates a nearly uniform 
1.4 T magnetic field in the z direction and provides pre- 
cision tracking and momentum measurement of charged 
particles within \r]\ < 1. The silicon tracker, which covers 
the 1 77 1 < 2 region, is composed of the innermost detector 
(LOO) 0, the Silicon Vertex Detector (SVXII) @, and 
the Intermediate Silicon Layers (ISL) [Hi] • LOO is a layer 
of single-sided radiation-hardened silicon strips mounted 
directly on the beam pipe at a radius ranging from 1.35 
cm to 1.62 cm. SVXII is an approximately 95 cm long 
cylinder of five layers of double-sided silicon microstrips 
covering a radial region between 2.5 cm and 10.7 cm. The 
ISL employs the same sensors as SVXII and covers the 
radial region between 20 cm and 28 cm, with one layer 
in the central region and two layers at larger angles. The 
COT [111 ], which spans 310 cm in length at a radial dis- 
tance ranging between 43 and 132 cm, contains four axial 
and four ±2° stereo superlayers of azimuthal drift cells. 
Axial and stereo superlayers alternate radially with one 
another. The COT provides full coverage in the |ry| < 1 
region, with reduced coverage in the region 1 < \r}\ < 2. 

Sampling calorimeters, divided into an inner elec- 
tromagnetic and an outer hadronic compartment, sur- 
round the solenoid. Except for limited areas of non- 
instrumented regions ( "cracks" ) , the calorimeters provide 
full azimuthal coverage within \rj\ < 3.6. All calorimeters 
are split into towers with projective geometry pointing at 
the nominal interaction vertex Q ■ Embedded in the elec- 
tromagnetic compartment, a shower maximum detector 
provides good position measurements of the electromag- 
netic showers and is used in electron identification [l2j]. 

The muon detection system consists of stacks of drift 
chamber modules backed by plastic scintillator counters. 
The stacks are four layer deep with laterally staggered 
cells from layer to layer to compensate for cell-edge in- 




FIG. 1: Elevation view of half of the CDF II detector, showing 
the inner microstrip detector, the Central Outer Tracker drift 
chamber, the electromagnetic and hadronic calorimeters, the 
muon drift chambers and scintillation counters. 



efficiencies. Four separate systems are used to detect 
muons in the \rj\ < 1.5 region. The central muon detec- 
tor (CMU) [13| is located behind the central hadronic 
calorimeter at a radius of ~ 3.5 m from the beam axis, 
covering the \rj\ < 0.63 region. The central muon upgrade 
detector (CMP) is arranged to enclose the \rj\ < 0.54 re- 
gion in an approximate four-sided box. It is separated 
from the CMU by the additional shielding provided by 60 
cm of steel. The central muon extension (CMX) extends 
the muon identification to the region 0.6 < \rj\ < 1.0. The 
more forward region (1.0 < \rj\ < 1.5) is covered by the 
intermediate muon detector (IMU). Table [T] summarizes 
the characteristics of the CDF sub-detectors used in this 
analysis. 

CDF uses a three-level trigger system to select events 
to be recorded on tape, filtering the interactions from a 
1.7 MHz average bunch crossing rate to an output of 75- 
100 Hz. This analysis uses data from triggers based on 
leptons with high transverse momentum Pt, as expected 
from the leptonically decaying W's in the event. The first 
two trigger levels perform limited reconstruction using 
dedicated hardware, which reconstructs tracks from the 
COT in the r — <f> plane with a transverse momentum 
resolution better than 2% x [GeV/c] fig . The electron 
trigger requires a coincidence of a COT track with an 
electromagnetic cluster in the central calorimeter, while 
the muon trigger requires that a COT track points toward 
a set of hits in the muon chambers. The third level is a 
software trigger which runs offline algorithms optimized 
for speed. 
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TABLE I: CDF II sub-detectors, purposes, resolutions or acceptances. 



Component 



Purpose 



Resolution / Acceptance 



Reference 



Silicon System 



Hit position 



Impact parameter 
Interaction Vertex Position 



11 nm (LOO) 
9 iim (SVXII) 
16-^23 (im (ISL) 
40 \im 
70 jim 



\7] 



COT 



Hit position 
Momentum measurement 



140 fim 
^ = 0.15 x P T [GeV/c] 



HI] 



Central Calorimeters 

Electro-magnetic calorimeter 
Shower Max Detector 
Hadron Calorimeter 
Wall Hadron Calorimeter 



Energy 
Position 
Energy 
Energy 



EM. 

E 



<Ze_ 

E 
"E 

E 



13.5%/ V E T [GeV] ffi 2% 
2 mm 



50.0%/ t /eJgW] ® 3% 
75.0%/ y/E[GeV] 4% 



[14] 
[15] 
[16] 

m 



Forward Calorimeters 

Electro-magnetic calorimeter 
Shower Max Detector 
Hadron Calorimeter 



Energy 
Position 
Energy 



"E 
E 



E 



16.0%/ yE[GeV] e 1% 

1 mm 
80.0%/ v 'E[GeV](B5% 



[17] 
[12] 

[7] 



Muon Systems 

CMU 
CMP 
CMX 



Muon Detection 



P T >IA GeV/c 
P T >2.2 GeV/c 
P T >l-4 GeV/c 



[13], [7] 

® 
[7] 



III. DATA SAMPLE 

The signature of tt dilepton events consists of two 
large transverse momentum leptons (e or fx), large miss- 
ing transverse energy ($r), two jets originating from b 
quarks, and possible additional jets from initial and final 
state radiation. We select dilepton events from inclusive 
high-.Pr electron and muon triggers using the standard 
CDF lepton+track algorithm, as described in the next 
sections. 

The main expected background processes in the dilep- 
ton sample are VF+jets with a jet misidcntified as a lep- 
ton ("fakes"), Drell-Yan (Z/j* -> e+e~, /x + M~, t + t~), 
and diboson events (WW, WZ, ZZ) with additional 
jets. In the case of Drell-Yan, non-physical $t can be 
faked by mis-measured jets or leptons. The contribution 
of these processes to the selected data sample is reduced 
by optimized selection cuts. 



A. Trigger 

A high-transverse- momentum lepton is required by the 
trigger. For a central electron candidate, an electro- 
magnetic calorimeter cluster with Et = E ■ sin 9 > 
18 GeV, accompanied by a matched COT track with 
Pt = P ■ sin# > 9 GeV/c, is required. For an elec- 
tron in the plug region (1.1 < < 2.0), the trigger re- 
quires an electromagnetic cluster in the calorimeter with 
E T > 20 GeV and E T > 15 GeV. For muon candidates 
two or more hits in the outer muon chambers matching 
a track of Pt > 18 GeV/c in the central tracker are re- 



quired. 



B. Leptons 

The LTRK selection aims at selecting two charged lep- 
tons of opposite charge with a greater acceptance than if 
tight lepton selection cuts were applied on both leptons. 
One lepton ("tight lepton") must have a well- measured 
track reconstructed from the interaction point with asso- 
ciated hits in the COT and SVX. For muons, the track is 
required to be compatible with hits in the muon chambers 
and to have Pt > 20 GeV/c and \r]\ < 1. For forward 
electrons a calorimetry-seeded tracking algorithm is used 
to identify tracks since the plug region is not well covered 
by the COT. In the case of electrons, the track is required 
to point to an electromagnetic cluster with Et > 20 GeV 
and < 2. Tight leptons must also satisfy an isolation 
requirement, i.e. the additional Et in a cone of radius 
AR = \J Ar] 2 + Atp 2 = 0.4 about the lepton trajectory 
must not exceed 10% of the lepton Et- 

The other lepton ("track lepton") is required to be a 
well-measured track originating at the interaction point 
with \r}\ < 1 and P T > 20 GeV/c. The track lepton 
must be isolated, which means that the ratio between the 
additional transverse momentum of tracks in a AR = 0.4 
cone around the track lepton and the overall Pt in the 
cone is less than 10%. Compared to the dilepton selection 
("DIL" [lj|) LTRK relaxes the calorimeter constraints on 
the track lepton in order to recover those events in which 
a lepton hits a detector crack. We refer to @ for a more 
detailed comparison between LTRK and DIL. 
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C. Jets 

Jets are the final products of quark hadronization. 
They are identified by looking for clusters of energy in the 
calorimeter using the JETCLU cone algorithm [2(|. The 
jet search is seeded by towers with Et > 1 GeV. Starting 
from the most energetic seed, all seeds within a 7x7 bins 
wide area around the seed are grouped into a cluster and 
the centroid is calculated. Seeds cannot belong to more 
than one cluster. All towers with Et > 0.1 GeV within 
a AR — 0.4 cone about the cluster centroid are added to 
the cluster and the centroid is recalculated. The proce- 
dure is iterated and a final step of splitting and merging 
is performed in order not to include the same tower in 
more than one jet. 

Jet transverse energy is corrected for non-uniformities 
in the calorimeter response as a function of jet rj, multiple 
pp interactions, and the hadronic jet energy scale of the 
calorimeter pl| . Events are required to have two or more 
jets with E T > 20 GeV and |r?| < 2. 

D. Missing Transverse Energy 

The definition of the uncorrected missing transverse 
energy is: 

% = -^E&ii (1) 

i 

where the sum is performed over all towers with a de- 
posited energy of at least 0.1 GeV. hi is the transverse 
unit vector pointing from the CDF geometrical center to 
th e I t h tower. 

IfiT is corrected to compensate for the following effects: 

• the interaction vertex displacement with respect to 
the CDF geometrical center: $t is recalculated 
with hi (Eq. [1]) as having the origin in the interac- 
tion point. 

• potential jet mis-measurement: if a track within 
the jet cone has a transverse momentum larger than 
the jet trasnverse energy, the difference between the 
Pt of the highest-Pr track and the jet Et is added 

to #r. 

• muons: to correct J^t for the identified muons and 
to account for their minimum ionization contri- 
bution in the calorimeters, the difference between 
muon calorimeter Et and muon Pt is added to $t- 

• jet corrections: J£t is updated according to the cor- 
rections applied to the jet transverse energies, as 
explained above. 

After corrections are applied the magnitude of the 
missing transverse energy is required to be larger than 
25 GeV. 



E. Final Selection Cuts 

Several topological vetoes are implemented in order 
to reduce the impact of backgrounds in the sample. 
Background contributions from Z boson decays yield- 
ing overestimated T^t are removed by raising the $t 
requirement to 40 GeV and the invariant mass of the 
tight lepton+track lepton pair to be inside the Z mass 
window ([76, 106] GeV/c 2 ). Large azimuthal separa- 
tions between the J^t and jets (A<fi > 25°), tight lepton 
(A0 > 5°), and track lepton (A0 > 5°, Ac/) < 175°) are 
required. These requirements have been implemented in 
order to reduce the number of events where mis-measured 
leptons or jets lead to overestimated Et, mostly con- 
tributed by the Drell-Yan process. A lower cut on the 
angle between the tight lepton and Et is applied to re- 
duce the acceptance for Z/j* — > p,fi as electron+track, 
where high-Pr muons are misidentified as electrons be- 
cause of the emission of bremsstrahlung photons. The 
requirement of a minimum azimuthal angle between jets 
and Et is dropped if Et > 50 GeV, since such large val- 
ues of missing transverse energy are not expected to arise 
from jet mis-measurements. 

Events with muons from cosmic rays or electrons origi- 
nating from the conversion of photons are removed. Cos- 
mic muons are identified by requiring a delayed coinci- 
dence of the particle hits in the calorimeter 22]. Con- 
versions are identified by pairing the electron track to 
an opposite sign track originating from a common ver- 
tex (23]. 



F. Sample Composition 

Table |TT] summarizes the ti signal and background rates 
expected for a LTRK sample corresponding to an inte- 
grated luminosity of 2.9 fb -1 . Depending on the pro- 
cess, background rates are estimated using simulated 
or data events. Simulated events are generated with 
the pythia |23l) Monte Carlo program, which employs 
CTEQ5L [2 4j parton distribution functions, leading- 
order QCD matrix elements for the hard process sim- 
ulation, and parton showering to simulate fragmentation 
and gluon radiation. A full simulation of the CDF II 
detector [25[ is applied. Diboson and Z/-f* — > t+t~ 
rates are estimated with simulated events, while Z/j* — > 
e + e~,[i + fi~ rates are estimated with a mixture of data 
and simulation. We use Z/j* — > e + e~,fi + fi~ simulated 
events to predict the ratio of events in different kine- 
matic regions, while we use data to normalize the over- 
all rates. The expected fakes from IF+jets and ti single 
lepton events with a jet misidentified as a lepton are esti- 
mated with VF+jets data [26[ . Signal acceptance and ex- 
pected rate are evaluated using simulated it events with 
a cross-section of 6.7 pb [13] and a top quark mass of 175 
GeV/c 2 . 
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TABLE II: Expected numbers of it signal and background 
events with statistical uncertainties for the LTRK data sam- 
ple. A it cross-section of 6.7 pb at a top quark mass of 
175 GeV/c 2 is assumed. 



Process 


Expected number 


Signal (tt) 


162.6 ± 5.1 


WW 


10.5 ± 1.0 


WZ 


3.8 ± 0.3 


ZZ 


0.9 ± 0.1 


Z/j* -* e + e~ 


20.8 ± 6.0 


Zh* n+fi- 


9.1 ± 3.1 


Z/y* -> t + t~ 


19.6 ± 2.4 


Fakes 


80.2 ± 15.7 


Total background 


145.0 ± 17.3 



IV. MASS RECONSTRUCTION 

In this section we describe the procedure to reconstruct 
an event-by-event preferred top quark mass (m[ eco ). In 
the next sections we will explain how the m 
tion is used to extract the top quark mass 



reco\ 
t 

r t eco distribu- 



A. Kinematics in the Dilepton Channel 

To reconstruct the ti event one needs to get 4-momenta 
for six final state particles, 24 values in total. These final 
state particles are two leptons and two neutrinos from 
W's decays, as well as two jets originated from the top- 
decay b quarks. Out of the 24 final quantities, 16 (jet 

and lepton 4-momenta) are measured, two ($t compo- 
nents in the transverse plane) are obtained by assum- 
ing overall transverse momentum conservation, and five 
constraints are imposed on the involved particle masses 
(m^- = mw+ = where mw = 80.4 GeV/c 2 [3j, 

m t — mj, m v — m— = 0) . The event kinematics is there- 
fore under-constrained. One must assume that at least 
one more parameter is known in order to reconstruct the 
kinematics and solve for the top quark mass. 



B. Neutrino <f> Weighting Method 

The method implemented in this work for reconstruct- 
ing the top quark mass event by event is called "Neutrino 
</> Weighting Method" . This method was previously de- 
scribed in [6j . In order to constrain the kinematics a scan 
over the space of possibilities for the azimuthal angles of 
the neutrinos (4> Ul , <j> V2 ) is used. A top quark mass is re- 
constructed by minimizing a chi-squared function (x 2 ) in 
the dilepton it event hypothesis. The x 2 has two terms: 



X 



v 2 - 
Areso 



Xc 



(2) 



The first term takes into account the detector uncer- 
tainties, whereas the second one constrains the parame- 



ters to the known physical quantities within their uncer- 
tainties. The first term is as follows: 



Xreso 



P l T? 



1=1 



E 



2 

!£>(^t/(i#|i£)) 



(3) 



<JEu' 



With the use of the tilde (~) we specify the parameters 
of the minimization procedure, whereas variables without 
tilde represent the measured values. £?tj are the transfer 
functions between b quark and jets: they express the 
probability of measuring a jet transverse momentum Pj, 

from a b quark with transverse momentum P^,. We will 
comment on ^ t f in Section IIV CI The sum in the first 
term is over the two leptons in the event; the second sum 
loops over the two highest-i?T (leading) jets, which are 
assumed to originate from the b quarks (this assumption 
is true in about 70% of simulated ti events [6J|). 

The third sum in Eq. [3] runs over the transverse com- 
ponents of the unclustered energy {Efj, E^), which is 
defined as the sum of the energy vectors from the towers 
not associated with leptons or any leading jets. It also 
includes possible additional jets with Et > 8 GeV within 
\rj\<2. 

The uncertainties (cp T ) on the tight lepton Pt used 
for identified electrons (e) and muons (/x) are calculated 
as 0: 



pe 
T 

7 Pt 



0.135 2 



P|[GeV/c] 



0.02 2 



= 0.0011 • P£[GeV/c] 



T 



(4) 



(5) 



The track-lepton momentum uncertainty is calculated 
as for the muons, since momentum is measured in the 
tracker for both electrons and muons. Uncertainty for the 
transverse components of the unclustered energy, ueu , is 

defined as 0.4 



£unclr 



where £^ ncl is the scalar 



'[GeV] [2 

sum of the transverse energy excluding the two leptons 
and the two leading jets. 

The second term in Eq. [2l Xconstn constrains the pa- 
rameters of the minimization procedure through the in- 
variant masses of the lepton-neutrino and of the lepton- 
neutrino-leading-jet systems. This term is as follows: 



Xc 



2\n{<? BW {mt: i \ 



m W, Twin/)) 
TmJ) 

r™J) 

rht is the parameter giving the reconstructed top quark 
ass. 3^Bw{m inv \ m,T) 



2]n(& BW (m l ^\ 
2m(JW(m^f' i3 |mt, 



(6) 



indi- 



cates the relativistic Breit-Wigner distribution function, 
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which expresses the probability that an unstable particle 
of mass m and decay width T decays into a system of 
particles with invariant mass m inv . We use the PDG Q 
values for mw and T mw . For the top width we use the 
function 



r - ^ F 



m 3 (1 _^ f) 2 (1 + 2 ^ ) 



(7) 



according to Ref. [29[. This new formulation of the 
Xconstr term helps to decrease the statistical error of the 
top mass reconstruction by 20%. 

The longitudinal components of the neutrino momenta 
are free parameters of the minimization procedure, while 

the transverse components are related to $t and to the 



assumed 



as follows: 



pv\ — 
± x — 


T 


■ cos(<f> Vl ) 

sin(0„ 2 )-$ 


Ty -cos(tj>„ 2 ) 






sill^^ — 




pv\ — 

y ~ 


pv\ 
T 

Vt x 


• sin^J 

sin(0„ 2 )-^ 


Ty -cos(0„ 2 ) 






sin(0,y 2 — 




pvi 

r x — 


T 


• cos(^ 2 ) 


Ty -COs(4> vl ) 






sin(0 M1 — 


<t>U 2 ) 


pvi — 

y ~ 


T 


• sin(^ 2 ) 

sin(</) 1/1 )-| 


Ty -cos(0„ 1 ) 






sin(<f} l ^ 1 — 


<f>»2 ) 



cos(0 !/1 ) 



sm((i) ui ) 



cos(</>„ 2 ) 



(8) 



sin 



The minimization procedure described above must be 
performed for all the allowed values of <j) L , 1 , 0„ 2 in the 
(0,27r) x (0, 2ir) region. Based on simulation, we choose 
a <j) Vl , cf> U2 grid of 24x24 values as inputs for the min- 
imization procedure. In building the grid we avoid the 
singular points at <f> vi = (f> U2 + k ■ 7r, where k is inte- 
ger. For these points, which correspond to a configu- 
ration where the two neutrinos are collinear in trans- 
verse plane, the kinematics of the event cannot be re- 
constructed using Eqs. [3J- M Avoiding these points in 
our procedure does not effect the reconstruction of the 
top mass central value, but rather effects the width of 
the mass distribution per event. Note from Eq. [8] that 
performing the transformation — > <p v + ir leaves P£ 
and Py unchanged, but reverses the sign of P£- We ex- 
clude unphysical solutions (Pp < and/or P^ 2 < 0) and 
choose the solution which leads to positive transverse mo- 
menta for both neutrinos. This decreases the number of 
grid points to 12x12. At each point 8 solutions can ex- 
ist, because of the two-fold ambiguity in the longitudinal 
momentum for each neutrino and of the ambiguity on the 
lcpton-jet association. Therefore, for each event, we per- 
form 1152 minimizations, each of which returns a value of 
and Xijk (h j = !,-■■, 12; k = 1, . . . , 8). We define 
4 • ln(r mt ), which is obtained by using Eq. 



l ijk 



X 



XI 



[5] where S?bw is substituted with -, — 2 r— rr-r . 

and select the lowest x' 2 solution for each point of the 



(<j) Vl , <f> V2 ) grid, thereby reducing the number of obtained 
masses to 144. Each mass is next weighted according to 



~X?J 2 



y 12 ^y 12 . 



(9) 



A top quark mass distribution is built in order to identify 
the most probable value (MPV) for the event. Based 
on a result of the simulation, the following procedure 
for improving the performance of solution- weighting was 
implemented. Masses below a threshold of 30% the MPV 
bin content are discarded, and the remaining ones are 
averaged to compute the preferred top quark mass for 
the event. 



C. Transfer Functions 

Since jet energy corrections have been calibrated on 
samples dominated by light quarks and gluons, we need 
an additional correction for a better reconstruction of 
the energy of 6-quark jets. In Equation [3l we intro- 
duced the transfer functions <5^t/, which allow us to 
step back from jets to partons. These functions of 
jet rj and Pr are defined as the parameterization of 
£ = ^p^-i uark _ p^ et yp^ et distributions, built from a 
large sample of simulated it events. The 6-quark jets 
in the simulation are recognized using true MC infor- 
mation. Jets with axis within a R = 0.4 cone about the 
generated b quarks are used. The influence of 6-quark P T 
spectra on the £ distributions is minimized by choosing 
the weights inversely proportional to the probability den- 
sity of p}^i uark _ Also, this greatly reduces dependence 
of the transfer functions on rut- 
in order to parameterize the above distributions we 
found the following expression to be adequate: 

7776 - n ere-Ti , ( «" 



-0.5(5 



'+exp(- 



V2tt7 2 
77(1 - 76 ) 



-0.5(5 



(1 



27T75 
" 77) 



(10) 



-0.5( 



j— ra \ 2 



27T73 



The parameters 71 • • • 7s are derived from the fit. The 
distributions are built for three regions: \i]\ < 0.7, 
0.7 < \n\ < 1.3, and 1.3 < \r]\ < 2.0. 

Figure [2] shows the distributions and the transfer func- 
tions for a number of (\r]\,P^ e ) regions. 10 GeV/c wide 
Pr bins are used from 30 GeV/c to 190 GeV/c for \r/\ < 
0.7, from 30 GeV/c to 150 GeV/c for 0.7 < \rj\ < 1.3, 
and from 30 GeV/c to 110 GeV/c for 1.3 < \r)\ < 2.0. A 
single bin is used above and below these regions. 



V. TOP QUARK MASS DETERMINATION 

The selected data sample is a mixture of signal and 
background events. In order to extract the top quark 
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FIG. 2: Examples of the transfer functions of b quarks into jets used in the fit. These functions of jet r\ and Pt are defined as 
the parameterization of (P T ~ 9uar 
the parameterization with Eq. 1101 



the parameterization of (P,£ quar — P^ et )/ P^ et distributions. The points are from the simulated tt events. The curves show 



mass, the reconstructed top quark mass distribution 
in data is compared with probability density functions 
(p.d.f.'s) for signal and background by means of a like- 
lihood minimization. P.d.f.'s are defined as the parame- 
terizations of m r t eco templates obtained by applying the 
neutrino <f> weighting method on simulated signal and 
background events, which are selected according to the 
lcpton+track algorithm. 



two Gaussian distribution functions, as: 



r> I reco\ s C lP6 -0.5( " f - P1 

P s \m t \m t ) = r— — e p 2 



+cxp(- 

- c. 

2np 2 

+ ci(l-Pe) e ~o.5( m '"°-^ ) 2 

( 1 _ C l) e -0.5( "'^°- C2 ) 2 
27TP3 



(11) 



A. Templates 

Signal templates are built from it samples generated 
with PYTHIA for top quark masses in the range 155 to 
195 GeV/c 2 in 2 GeV/c 2 steps. They are parameterized 
in a global fit by using a combination of one Landau and 



P s , the signal p.d.f., expresses the probability that a mass 
m reco j g rec0 nstructed from an event with true top quark 
mass nit- The constants c\ and ci are set a-priori to 
adhere to the features of the template shape. The pa- 
rameters pi , . . . , Pq depend on the true top quark mass 
m t are calculated as: 

Pk = a k + a k+6 -(m t [GeV/c 2 ]-l75) fc=l,...,6 (12) 
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The parameters ak are obtained from the fit to the signal 
templates. Figure shows a subset of templates along 
with their parameterizations (solid lines). 

A representative background template is built by 
adding fakes, Drell-Yan, and diboson templates. These 
templates have been normalized to the expected rates 
reported in Table |TTJ The fakes template is built from 
W^+jets data events by weighting each event according 
to the probability for a jet to be mis-identified as a lep- 
ton (fake rate) [2y]. Drell-Yan and diboson templates 
are built from samples simulated with PYTHIA and ALP- 
GEN [3(| + PYTHIA respectively. The combined back- 
ground template is fitted with a sum of two Landau and 
one Gaussian distribution functions, as: 



and 



~3i >> 



fcl(l ~ Ai) c -0.5( m '"°"^ 4 ) 2 



2tt(3 5 



(13) 



(1 - fcl) o -0.5( m '"° fc2 +cxp(-^ 



TT 



2tt/3 3 



where the fit parameters 0i ■ ■ ■ fa are m t -independent. 
The constants k\ and hi are set a-priori to adhere to 
the features of the template shape. The combined back- 
ground template and its parameterization (solid line), 
Drell-Yan, diboson, and fakes templates are plotted in 
Figure |U 



Jfparam = exp{-0.5[(d- a ) T U 1 (a - do) 
+ 0-fo) T V- 1 0-&)]}. 



(17) 



The top quark mass estimator (mj op ) returned by the 
minimization is the mass corresponding to [— lnJ§f] m j n . 
The shape likelihood term, J£ s ha P e (Eq. H5]l . expresses 
the probability of an event being signal with the top 
mass mtop or background. The signal (P s ) and back- 
ground (P(0 probabilities are weighted according to the 
number of signal (n s ) and background (rib) events, which 
are floated in the likelihood fit. In the fitting procedure, 
rib is constrained to be Gaussian-distributed with mean 
value n e b xp and standard deviation u n ^v , as shown in Eq. 
I16[ while (n s + rib) is the mean of a Poisson distribu- 
tion of N selected events. In this manner, the number 
of signal events is independent of the expected tt lep- 
ton+track events in a particular assumption of the tt 



cross-section value. J£ m 



constrains the parameters 



of the signal (a) (see Eq. [T2)) and background (/?) (see 
Eq. [T3|) p.d.f.'s. These p.d.f.'s have Gaussian distribu- 
tion with mean values (do) and Q3o) obtained from the 
signal and background templates fit. U and V are the 
corresponding covariant matrices for a and (3 returned 
from the MINUET [3l| minimization. 



VI. CALIBRATION OF THE METHOD 



B. Likelihood Minimization 

The top quark mass estimator is extracted from the 
data sample by performing an unbinncd likelihood fit 
and minimization. The likelihood function expresses the 
probability that a m r t eco distribution from data is de- 
scribed by a mixture of background events and dilcpton 
tt events with an assumed top quark mass. Inputs for the 
likelihood fit are the reconstructed mass (m n ), the sim- 
ulated signal and background p.d.f.'s, and the expected 
background. The background expectation (nf xp =145.0) 
and its uncertainty (a n ^p = 17.3) are taken from Table [TT1 
The likelihood takes the form: 



^shape ' -^backgr ' -^param^ 



(14) 



where 



shape 



_ e -(«=+«<.) . (n s +n b ) N 
~ AH 

N v ( \ \ I u ( \ ( 15 ) 

tt n s ■ P s \m n \m t0 p) + n b ■ Pb(m n ) 

, n s + n b 



backgr 



exp(- 



-{rib 



exp\2 



2a 2 



(16) 



The method described above is calibrated in order to 
avoid systematic biases in the measured top quark mass 
and in its uncertainty. Calibrations are performed by 
running a large number (10 4 ) of "pseudo-experiments" 
(PE's) on simulated background and signal events where 
the true top quark mass is known. Each PE consists of 
determining the number of signal [Nf ) and background 
(Nf ) events in the sample, drawing N^ E masses from 
a signal template and N b E from the background tem- 
plate, and fitting the mass distribution to a combination 
of signal and background p.d.f.'s, as described in Section 
IVl A top quark mass (m{ lt ) and its positive and neg- 
ative statistical uncertainties (a + and cr~ ) are returned 
by the fit. Numbers of signal and background events are 
generated according to Poisson distributions with means 
given in Table HT1 

For each input top quark mass the median of the m t %t 
distribution is chosen as the mass estimate (m° ut ). The 
distributions of m° versus input mass (to*) and the bias, 
defined as AM = m° ut — m t , are shown in Figure [3 The 
uncertainty bars are determined by the limited statistics 
of the signal and background templates. Both fits in Fig- 
ure[S]are performed in the mass range 159-191 GeV/c 2 . 
The slope of the straight line in the upper plot is consis- 
tent with one, while the average bias (horizontal line in 
the lower plot) is -0.13±0.10 GeV/c 2 . Although this 
value can be considered compatible with zero within un- 
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159GeV/c 2 165GeV/c 2 171 GeV/c 2 





FIG. 3: Signal templates and fitting functions (solid lines) for a number of generated top quark masses. The parameterization 
is defined in Eq. 1 1 1 1 



certainties, we apply a shift of +0.13 GeV/c 2 to the result 
on data. 

In order to check the bias on the statistical uncertainty 
we use pull distributions, defined as follows: 



Jit 



pull 



nit 



(18) 



where a' = { 



The positive and negative 



if < mt 

\(T~ | if m( rt > nit 

statistical uncertainties are returned by MINUIT (rou- 
tine MINOS) [3l[. For each generated top quark mass, 
pull distributions are fitted by Gaussian functions (some 
examples are shown in Figure [5]) ■ 

The mean and width of the pull distributions versus 
generated top quark mass are shown in Figure [7J Error 
bars account for the limited statistics of signal and back- 
ground templates. The average width of pull distribu- 
tions is 1.009 ± 0.005. A width larger than one indicates 
an underestimate of the statistical uncertainty. Accord- 
ingly, the statistical uncertainty obtained from data is 
increased by a factor 1.009. 



VII. RESULTS 

The data sample used in this measurement corresponds 
to an integrated luminosity of 2.9 fb^ 1 . A total of 328 
LTRK candidates are found in data. Selected events are 
reconstructed and an experimental mass distribution is 
built. The likelihood constrained fit described in Section 
IV Bl is performed and the following estimate of the top 
quark mass with statistical uncertainties is obtained: 



m top 



165.35±|;|| GeV/c 



(19) 



The experimental top quark mass distribution is shown in 
Figure [HI The constrained fit returns 181.4^2X3 signal 
and 146.1^5 J background events. The observed rates 
are in good agreement with expectations (Table HTf . 

As a check, we remove the Gaussian constraint on the 
number of background events in Formula 1 141 The uncon- 
strained fit returns: 



m t c 



165.33 



+3.39 
^3.28 



GeV/c 2 



(20) 



with 178.6_3^i signal and 149. background events. 
The top quark mass and the number of signal and back- 
ground events from unconstrained and constrained fits 
are in agreement. 
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FIG. 4: Drell-Yan, fakes, diboson and combined background templates. The fitting function (solid line), defined in Eq. I13I is 
superimposed to the combined template. 
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FIG. 5: Results from pseudo-experiments. The upper plot 
shows m° ut versus input masses, while the lower one shows 
the bias. 
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FIG. 6: Results from pseudo-experiments: pull distributions 
for generated mass samples at rat = 175 GeV/c 2 (left) and 
m t = 181 GeV/c 2 (right). Distributions are fitted to Gaus- 
sian functions (solid line), returning the indicated means and 
standard deviations. 



The top quark mass and its statistical uncertainty ob- 
tained from the constrained fit (Eq. IT5|) are corrected for 
the expected systematic 0.13 GeV/c 2 shift, and for the 
1.009 width of the pull distribution (Section IVI[) respec- 
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FIG. 7: Results from pseudo-experiments: mean and width 
of the pull distributions versus generated top quark mass are 
shown in the upper and lower plots respectively. 
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FIG. 9: Expected statistical uncertainties from pseudo- 
experiments generated with a top quark mass of 165 GeV/c 2 . 
The arrows indicate the uncertainties found in this measure- 
ment. 
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tively. The final value is: 

m top = 165.5t|;^(stat.) GeV/c 2 



(21) 



100 150 200 250 300 350 
Reconstructed mass (GeV/c 2 ) 

FIG. 8: Two-component constrained fit to the 328-event 
TTRK data sample. Background (dark gray) and sig- 
nal+background (light gray) p.d.f.'s, normalized according to 
the numbers returned by the fit, are superimposed to the 
reconstructed mass distribution from data (histogram). The 
insert shows the fitted mass-dependent negative log-likelihood 
function. 



In order to check that the measured statistical uncer- 
tainty is reasonable, a set of PE's is performed on 
simulated background and signal events with m t = 
165 GeV/c 2 (close to the central value of the constrained 
fit), as explained in Section [VTl The obtained positive 
and negative error distributions along with the observed 
values (arrows) are shown in Figure O We found that 
the probability for obtaining a precision better than that 
found in this experiment is 82%. 



VIII. SYSTEMATIC UNCERTAINTIES 

Since our method compares findings to expectations 
estimated from Monte Carlo simulations, uncertainties 
in the models used to generate events cause systematic 
uncertainties. Other systematic uncertainties arise from 
the potential mis-modeling of the background template 
shape. 

The procedure for estimating a systematic uncertainty 
is as follows. The parameters used for the generation of 
events are modified by ± 1 standard deviation in their 
uncertainties and new templates are built. PE's from the 
modified templates are performed using the same p.d.f.'s 
as in the analysis. The obtained medians of the top quark 
mass distribution from PE's and the nominal top quark 
mass are used to estimate the systematic uncertainty. 
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The source of each systematic uncertainty is assumed to 
be uncorrelated to the other ones, so that the overall sys- 
tematic uncertainty is obtained by adding in quadrature 
the individual uncertainties. The systematic uncertain- 
ties along with the total uncertainty are summarized in 
Table fVl In the following, we describe how each system- 
atic uncertainty is evaluated. 



A. Jet energy scale 

The measured jet energy is corrected according to 
the measured and simulated calorimeter response to 
electrons and hadrons [32j. Jet corrections also correct 
for the non-uniformities in calorimeter response as a 
function of \rj\, effects of multiple pp collisions, the 
hadronic jet energy scale, deposited energy within the 
jet cone by the underlying events, and out-of-cone jet 
energy lost in the clustering procedure. The systematic 
uncertainty due to the jet energy scale ( JES) is estimated 
from signal and background events in which each jet 
energy correction has been shifted by ± 1 standard 
deviation in the energy scale factor. Shifted signal and 
background templates are built and two sets of 10 4 PE's 
are performed. The systematic uncertainty for each level 
of corrections is taken as (mf — m t ~)/2, where mf and 
m~[ are the top quark masses found respectively for a 
lower and upper shift of the parameter. The individual 
uncertainties are summed in quadrature in order to 
obtain the JES systematic uncertainty. Results are 
reported in Table IIIII The systematic uncertainty in the 

TABLE III: Summary of jet energy scale systematic uncer- 
tainties 



Source 


Uncertainty (GeV/c 2 ) 


77 calorimeter non-uniformity 


0.6 


multiple interactions 


0.0 


hadronic jet energy scale 


2.2 


underlying event 


0.2 


out-of-cone energy loss 


1.8 


Total 


2.9 



top quark mass due to the JES uncertainty is 2.9 GeV/c 2 . 

Since jet energy corrections are estimated with stud- 
ies dominated by light quarks and gluon jets, additional 
uncertainty occurs on the 6-jet energy scale because of 
three main reasons [28| : 

1. uncertainty in the heavy- flavor fragmentation 
model; 

2. uncertainty in the fe-jet semileptonic branching ra- 
tio; 

3. uncertainty in the calorimeter response to energy 
released by 6-jets. 



The effect of the fragmentation model on the top quark 
mass is evaluated by reweighting events according to two 
different fragmentation models from fits on LEP [33| and 
SLD [3i| data, while effects of the uncertainties on the 
semileptonic fr-jet branching ratio (BR) and 6-jet energy 
calorimeter response are estimated by shifting the BR 
and the 6-jct energy scale. In all cases shifted templates 
are built and PE's are performed. The resulting shifted 
masses are used to estimate the systematic uncertainty 
due to each of the sources. These uncertainties are added 
in quadrature. The total systematic uncertainty in the 
&-jet energy scale is 0.4 GeV/c 2 . 

B. Lepton energy scale 

The uncertainty on the lepton energy scale may affect 
the top quark mass measurement. This uncertainty is 
studied by applying a ±1% shift to the Pt of leptons [2lj . 
Shifted templates are built and PE's are performed. Half 
of the difference of the resulting masses is taken as the 
systematic uncertainty on the top quark mass due to the 
lepton energy scale uncertainty. The systematic uncer- 
tainty in the lepton energy scale is 0.3 GeV/c 2 . 

C. Monte Carlo event generation 

Several systematic uncertainties are due to the model- 
ing of tt signal events. 

1. Monte Carlo generators 

The effect of the choice of a particular Monte Carlo 
generator is studied by comparing our default pythia 
generator to herwig. These generators differ in the 
hadronization models, handling of the underlying pp 
events and of the multiple pp collisions in the same bunch 
crossing, and in the spin correlations in the production 
and decay of tt pairs (implemented in HERWIG only) [35[ . 
The difference between masses obtained from sets of PE's 
performed with the two generators is found. The sys- 
tematic uncertainty due to our choice of Monte Carlo 
generators is 0.2 GeV/c 2 . 

2. Initial and final state radiation 

The effect of the initial and final state radiation (ISR 
and FSR) parameterization is studied, since jets radiated 
by interacting partons can be misidentified as leading jets 
and affect the top quark mass measurement. The system- 
atic uncertainty associated with ISR is obtained by ad- 
justing the QCD parameters in the DGLAP [36| parton 
shower evolution in tt events. The size of this adjustment 
has been obtained from comparisons between Drell-Yan 
data and simulated events [281 ]. Since the physical laws 
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that rule ISR and FSR are the same, the parameters that 
control ISR and FSR are varied together (IFSR) . Half of 
the difference in top quark mass from PE's performed 
on samples with increased and decreased IFSR is taken 
as the systematic uncertainty for the radiation modeling. 
The systematic uncertainty due to uncertainties in the 
initial and final state radiation is 0.2 GeV/c 2 . 



3. PDFs 

The uncertainty in reconstructing the top quark mass 
due to the use of sets of parton distribution function 
(PDF) comes from three sources: PDF choice, PDF pa- 
rameterization, and QCD scale {Aqcd)- The uncertainty 
due to the PDF choice is estimated as the difference be- 
tween the top quark mass extracted by using CTEQ5L 
(default) and MRST72 [jtfi]. The uncertainty due to 
PDF parameterization is estimated by shifting by ±1 
standard deviation one at a time the 20 eigenvectors of 
CTEQ6M [H]. Half of the differences between the shifted 
masses derived from PE's are added in quadrature. The 
measured mass differences between MRST72, generated 
with Aqcd = 300 MeV, and MRST75, generated with 
Aqcd = 228 MeV, H3] is taken as the uncertainty due to 
the choice of Aqcd- These systematic uncertainties are 
added in quadrature. Results are summarized in Table 
TVl The total systematic uncertainty due to uncertainties 
in the PDFs is 0.3 GeV/c 2 . 
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FIG. 10: Results from pseudo-experiments performed using 
events selected according to the number of interactions. 



to derive the systematic uncertainty. We multiply 0.26 



GeV/c /interaction by < 



data 



< N: 



data 



>=2.07 and < N™9 



> - < N*?c >, where 
>=1.50 are the average 



number of vertices in the selected data sample and simu- 
lated sample respectively. We obtain a 0.15 GeV/c 2 top 
mass uncertainty due to the event pile-up. 



D. Background template shape 



TABLE IV: PDF systematic uncertainties on top quark mass. 
The total systematic uncertainty is the sum in quadrature of 
the individual contributions. 



Source 


Uncertainty (GeV/c A ) 


PDF parameterization 


0.2 


PDF choice 


0.1 




0.2 


Total 


0.3 



4- Luminosity profile (event pileup) 

Pseudo-experiment simulations have only been made 
for a probability of multiple interactions in a single bunch 
crossing as appropriate for the collider luminosity during 
the first period of data taking (1.2 fb _1 integrated lumi- 
nosity.) A possible discrepancy between simulation and 
data collected at later times at higher luminosity may af- 
fect the top quark mass measurement. We evaluate this 
effect by running batches of PE's on ti events, selected 
according to the number of interaction vertices found in 
the event. 

The results from PE's are plotted against the num- 
ber of interactions and a linear fit is applied (Figure [TO")) . 
Since we do not see a significant mass dependence, we use 
the uncertainty (0.26 GeV/c 2 /interaction) on the slope 



The systematic uncertainties due to the potential mis- 
modeling of the background template shape were also 
estimated. We identify three independent sources for 
this systematic uncertainty: background composition, 
W+jets fakes shape, and Drell-Yan shape. The effect 
of the diboson shape is neglected because of the small 
expected rate of this background (Table QT]) . 

In order to estimate the systematic uncertainty for the 
background composition, fakes, diboson, and Drell-Yan, 
the expected rates are alternatively varied by plus or mi- 
nus one standard deviation (Table [H]) without changing 
the total number of expected background events. Half of 
the differences between ± 1 a shifted masses derived from 
PE's are added in quadrature. The systematic uncer- 
tainty due to uncertainty in the background composition 
is 0.5 GeV/c 2 . 

The uncertainty on the shape of the fake background 
template (Section IV A[) is modeled. The fake rate Et 
dependence is varied according to the fake rate uncer- 
tainties in each Et bin. Two shifted background tem- 
plates are built and used for PE's. The corresponding 
shift in mass is taken as the systematic uncertainty due 
to potential mis-modeling of fake shape. The top mass 
uncertainty due to uncertainty in the fake shape is 0.4 
GeV/c 2 . 

Drell-Yan events with associated jets can pass the se- 
lection because jet mis-measurements can cause a large 
unphysical $t- Mis- modeling of this effect is studied, 
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since it may affect the top quark mass measurement. Two 
modified Drcll-Yan templates are built by re-weighting 
Z/7* — > e + e~,/i + /i _ events. The weight has been opti- 
mized by looking at discrepancies in $t between Monte 
Carlo simulation and data. Results of PE's performed 
with the modified Drell-Yan templates are used to esti- 
mate the systematic uncertainty due to the possible fluc- 
tuation in the shape of this background. The mass sys- 
tematic uncertainty due to uncertainties in the shape of 
the Drell-Yan background is 0.3 GeV/c 2 . 



TABLE V: Summary of systematic uncertainties on the top 
quark measurement. 



Source 


Uncertainty 


Jet energy scale 


2.9 


fe-jet energy scale 


0.4 


Lepton energy scale 


0.3 


Monte Carlo generators 


0.2 


Initial and final state radiation 


0.2 


Parton distribution functions 


0.3 


Luminosity profile (pileup) 


0.2 


Background composition 


0.5 


Fakes shape 


0.4 


Drell-Yan shape 


0.3 



(GeV/c") 
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IX. CONCLUSIONS 

Using the template technique on a lepton+track sam- 
ple we measure a top quark mass of 

m top = 165.5^ 3(stat.)±3.1(syst.) GeV/c 2 

or (22) 
m top = 165.5t^GeV/c 2 

This result agrees with the world average top quark 
mass (mtop — 172. 4± 1.2 GeV/c 2 |38(), obtained by com- 
bining the main CDF and D0 Run I (1992-1996) and Run 
II (2001-present) results. 

Compared with our previous result (mt op — 169.7±9.8 
GeV/c 2 @ ), obtained on a / Ldt = 340 pb" 1 data sam- 
ple, a significant improvement in the total uncertainty 
has been achieved. The improvement due to the novel- 
ties in the analysis technique is estimated from PE's to be 
about 20%. The improvements which made this progress 
possible are the introduction of relativistic Breit-Wigner 
distribution functions in event reconstruction, along with 
TO top -dependent top width, while in [(J Gaussian distri- 
bution functions and a constant top width were used. A 
new feature of this analysis is the use of a larger statistics 
lepton+track sample which overlaps by only ~ 45% with 
the often used dilepton sample @ . 
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